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A description is given of three digital methods which have been used to recover ampli- 
tude and phase information from a modulated sinusoidal signal sampled at equal intervals 
of not more than one-sixth of a period. The first method, the "zero-crossing" met hod, is 
economical of computer time and, for modulation which is not too deep and does not contain 
frequencies near the carrier frequency, accurately recovers the phase and amplitude modu- 
lation. The second method, the "filter" method, is more laborious but it gives better accuracy 
and will operate with deeper and more rapid modulation. The third method, a statistical 
approach, will work with severly overmodulated signals, but it yields only a statistical sum- 
mary of the modulation. The methods were designed specifically for analysis of radio-star 
scintillation records but they may be applied to many other modulated signals. 



1. Introduction 

Three digital methods are described which have 
been used to recover amplitude and phase informa- 
tion from a modulated sinusoidal signal sampled at 
equal intervals of not more than one-sixth of a 
period. Although these methods were designed 
specifically for the purpose of analyzing the records 
described in the preceding paper, " Amplitude and 
angular scintillations of the radio source Cygnus-A 
observed at Boulder, Colorado", they may be ap- 
plied to many other modulated signals. In general, 
the desired modulation is imposed on a sinusoidal 
carrier of constant frequency and amplitude. The 
sample interval and the carrier frequency must be 
chosen carefully with regard to the expected range 
of modulating frequencies, in order to avoid biasing 
the results either by omission of or folding at fre- 
quencies of interest. 

All three of the methods described herein were 
designed for use with the IBM-650 digital computer 
and its peripheral equipment which included an 
automat c graph plotter. The input data were 
digital samples of the modulated carrier, sampled 
to an accuracy of \){ percent of full scale (six bits) 
and punched on paper tape. The first method of 
analysis, the " zero-crossing method", uses the prop- 
erties of the signal only at the times when it passes 
through its average value. Phase modulation is 
derived from variations in the intervals between 
these times; amplitude modulation is derived from 
the slope of the signal at these times. The second 
method, the "filter method", passes the signal 
through appropriate digital filters which recover the 
phase and amplitude modulation. The last method 
is a statistical one, used to analyze records which, 
because of overmodulation in amplitude and/or 
phase, are not amenable to treatment by either the 
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zero-crossing or filter methods. This statistical 
method cannot give amplitude and phase as a func- 
tion of time, but produces results which describe 
some statistical properties of these quantities. 

2. Zero-Crossing Method 

In this procedure the computer must deal with 
the modulated carrier at the times when it crosses 
its average value. The first step is to replace the 
original samples by deviations from their average 
value. Since drifts in the digital recording system 
might cause the carrier to undergo slow changes in 
average (d-c) level, changes unrelated to the modu- 
lation, it is usually better to use a moving average 
several periods in length than to use an average over 
the entire record. If a moving average is used, the 
process of replacing the original samples by devia- 
tions from the average is, in fact, a digital high-pass 
filter which removes the d-c level and its slow varia- 
tions, but does not appreciably affect the carrier or 
its modulation. For the radio-star observations, a 
moving average 30 samples in length (approximately 
5 periods) proved satisfactory. 

After the filtering operation, the computer can 
find the points of interest by searching for adjacent 
samples with opposite signs. Each such "crossing 
pair" of samples yields, by linear interpolation, a 
good estimate of the exact time of zero crossing. 
In order to extract phase information, first differ- 
ences are taken of the list of times of zero crossings. 
Each of these intervals is then replaced by its devia- 
tion from the mean interval. The cumulative sums 
of these deviations provide a list of numbers called 
the phase table. The phase table is simply a tabu- 
lation of instantaneous phase variations at half- 
period intervals. 

From the crossing pairs the computer derives 
another set of numbers called the amplitude table. 
This tabulation contains values proportional to the 
amplitude of the modulating signal. It is the sum 
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of the absolute values of the two samples involved 
in each crossing pair, and measures the instantaneous 
amplitude at half-period intervals. It will be 
demonstrated later that this simple method, which 
might be expected to be inaccurate, is, in fact, an 
effective means for extracting the envelope of the 
sinusoidal trace. Figure 1 gives examples of actual 
amplitude and phase plots, made by an automatic 
plotter from the amplitude and phase tables. Each 
plot represents 600 points (10 min) of data and has 
been expanded to occupy the full scale. The full- 
scale range of each phase plot is shown in electrical 
degrees. The reliability and sensitivity, together 
with a brief description of the drawbacks, of the 
zero-crossing method will be discussed in the follow- 
ing sections. 

3. Filter Method 

Amplitude and phase tables may also be obtained 
by passing the data through a system of digital 
filters. This analysis procedure is based upon a 
method suggested* by N. R. Goodman [1960], of 
Space Technology Laboratories. 

Consider a modulated sine wave which is sampled 
once every t seconds to produce a time sequence of 
n data points, which we designate x t . We are inter- 
ested in the amplitude (i.e., envelope), and phase 
of this record. A possible representation of the 
amplitude A(t) and phase a(t) is 

A(t) exp [ia(t)] = X(t)+iY(t) 

where X(t) and Y(i) are the real and imaginary parts 
of the complex representation of the power and 
phase. We may let the x t represent a set of sampled 
points from the X(t) process. If we can somehow 
construct from the x t a set of y t which will be a good 



representation of the Y(t) process, then we may 
compute A(t) and a(t) from the relations 

A(t) = [X 2 (t) + Y 2 (t)] l/ > 

a(t) = sin~ 1 [F(*)/il(*)]=cos- 1 [X(t)IA{i)l 

One possible representation of the y t in terms of 
the x t is 
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where the a s have suitable values. There are several 
ways in which one might evaluate the a s . Goodman 
shows that for the present purpose the a s should be 
chosen in such a way that, over the frequency range 
of interest, the spectral density matrix of the result- 
ing (x u y t ) process is a good approximation to the 
spectral density matrix of the actual [X(t), Y(t)] 
process. Following this line of thought, a set of a s 
were chosen to satisfy that criterion. For the case 
where the sampling rate was about 6 samples per 
period, the following a s proved satisfactory. 



do = 


0.16 


a\ = 


.15 
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OCTOBER 28,1958, CYGNUS-A, ELEVATION ANGLE = 18? DEPTH INDEX = 4,I08 Mc/s 





SEPTEMBER 21 , 1958 , CYGNUS- A , ELEVATION ANGLE =38°, DEPTH INDEX=2 108 Mc/s 
Figure 1. Examples of amplitude and phase plots obtained from the computer. 

Each record is 600 samples (10 min) in length. The bottom edge of the amplitude plot represents zero signal, and the plot has been normalized 
full scale. The scale of the phase plot, a, given on the figure, is in electrical degrees. 
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These a 8 were then used in conjunction with (1) to 
evaluate t he y t . 

Since the //, are a moving linear combination of the 
x u they are in fact the result of passing the x t through 
a linear digital filter. This filter (■?/ filter) has a 
narrow pass band, approximately unit gain, and 
7r/2 phase shift in its pass band. r \ ne phase shift, 
uniform over the pass hand, is accomplished by mak- 



ing the filter anti-symmetric; i.e., a s 



and 



a Q =0. Noise outside the pass band of the filter 
may be removed from the x t by means of a similar 
filter (x filter) with zero phase shift. This will be 
a symmetric filter, a s =a_ s . 

The amplitude and phase tables produced by either 
the "zero-crossing" or "filter" methods may be used 
for such subsequent purposes as spectral analysis 
and cross correlation between amplitude and phase. 

In order to test the ability of the "zero-crossing" 
and "filter" methods to reproduce accurately the 
amplitude and phase fluctuations inherent in the 
data, an artificial record with known amplitude and 
phase variations was analyzed using each method. 
The known amplitude and phase variations are 
shown by the heavy curves in figures 2 and 3. The 
light curves labeled "two-filter method" show the 
results of using both the x and y filters. To assess 
the possibility of noise contamination, the same 
artificial record was analyzed using only the y filter. 
The result of this analysis is shown in the figure 
labeled "one-filter method." 

Although both methods reproduce the known 
amplitude and phase variations quite well, the one- 
filter method includes a high-frequency component 
which is not evident in the two-filter method, indi- 
cating the existence of noise outside the pass band 
of the filter. 

The heavy curves (input) are tilted slightly with 
respect to the light curves (output ). This is because 
the input data is passed through a high-pass filter 
which removes trends. 

The figures labeled "zero-crossing method" show 
the corresponding results produced by the zero- 
crossing method when it is applied to the artificial 
data. 

The principal difference between the filter method 
and the zero-crossing method is that the former 
utilizes consecutive data points whereas the latter 
uses only pairs of consecutive points which have 
opposite signs (crossing pairs). For this reason the 
filter method has a twofold advantage: (1) mod 
2-7T ambiguities are less likely to occur between 
consecutive points than between crossing pairs, and 
(2) higher frequency components are preserved by 
the filter method since its uses data points which 
are spaced more closely in time. The two-filter 
method involves about four times as much computing 
time as does the zero-crossing method, but the 
fidelity of the more lengthy method is better. 

4. Statistical Method 

Both of the foregoing methods have been applied 
to artificial data containing various amounts of 
amplitude and phase modulation. Although the 



filter method is more capable of producing accurate 
amplitude and phase plots than is the zero-crossing 
method, both methods fail for some records. These 
highly modulated records contain significant phase 
and amplitude modulation with frequencies com- 
parable to the carrier frequency. 

In order to analyze such records, a statistical ap- 
proach was taken. This procedure does not give 
detailed amplitude and phase plots but instead gives 
two quantities, P n and a d , related respectively to 
the magnitudes of the amplitude and phase fluctua- 
tions. 

P n is evaluated by firs! finding the points on the 
actual modulated sine wave, P(f) } which lie above 
the peaks of a sine wave having the same average 
amplitude. The average value of these points, 
divided by the average amplitude, is defined as P n . 
The sampled points which enter into the evaluation 
of P n all result from amplitude modulation, since 
phase modulation can never increase the instantane- 
ous amplitude of the record above the peaks of the 
equivalent sine wave. Thus, P n is a statistical 
measure of the amplitude modulation on the record 
and is unaffected by the phase modulation. The 
standard deviation, a d , of the intervals between zero 
crossings, is a measure of the phase modulation alone 
since amplitude modulation cannot affect the time 
of zero crossings. 

It is desirable to relate the quantities P n and a d to 
quantities winch are more commonly associated with 
amplitude and phase effects. For example, the 
quantities 
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and a a were of interest in the preceding paper. 
Here A is the amplitude of the signal, A the mean 
amplitude, and <j a is the standard deviation of the 
phase fluctuations, a. 

One possible 4 way to relate P n and a d to (AA/A) 2 
and a a respectively, is to compute (AA/A) 2 and a a 
from the variations in amplitude and phase for a 
particular record as given by the "zero-crossing" or 
filter" method, and then to find P„ and a d for this 
same record. If enough records are treated in this 
manner, it is possible to obtain regression curves of 
(AA/A) 2 versus P n and of a a versus <r d . This pro- 
cedure not only gives the desired relations but yields 
confidence limits for P n and <j d as measures of (AA/A) 2 
and <r a respectively. However, it cannot be applied 
to the highly modulated records for which the sta- 
tistical quantities P n and <j d were originally intro- 
duced. A better method is to compute P n and <j d 
from artificial records which have a predetermined 
value of (AA/A) 2 and <r a , for the particular example 
under consideration . 

The artificial records were constructed in the fol- 
lowing way. A pair of random normal deviates (m u 
Ui) are used to compute the quantity 

P i =(a+m i ) 2 +n 2 i 

where a is a constant. P t may be thought of as the 
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square of the envelope of a signal at time t t which 
consists of the sum of a sine wave with amplitude a 
and noise with mean square amplitude of unity. 
The Pi are then passed through a low-pass filter 
which auto-correlates them over a length of 30 points. 
The artificial data points, d u are then computed 
from the following expression: 

di= P f sin (uoti+at), 



where a t , the phase fluctuation, is normally dis- 
tributed with mean zero, standard deviation a a and 
auto-correlation length of 30 ; and a> is some constant 
angular frequency, approximately 7r/3 per sample 
interval. The d t are then normalized so that the 
artificial records will have the same range in value 
of the d t for any a. 

The probability distribution function of the (P*) 172 
constructed in this manner has been derived by S. 0. 
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Figure 2. An illustration, using artificial data, of the effectiveness of the zero-crossing and filter methods in reproducing 

amplitude fluctuations inherent in the original record. 
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Rice [1944 and 1945]. (The reason for constructing 
the artificial data in this particular way was, of 
course, dictated by knowledge of the sort of signal to 
be analyzed. Signals with different modulation 
characteristics would require different techniques for 
producing the artificial data.) 



The method involving artificial data has the ad- 
vantage that it cannot be biased by mod 2w am- 
biguities since a a is a predetermined quantity. Thus 
the artificial records may be used to extend empiri- 
cally the relation between <r d and <r a to regions where 
<r a is much greater than 1 radian. P n) measured on 
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Figure 3. An illustration, using artificial data, of the effectiveness of the zero-crossing and filler methods in reproducing phase 

fluctuations inherent in the original record. 
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0.30 



Figure 4. The relationship, derived from artificial data, 
between the statistical quantity P n and the actual value of 

(Ai/i)2. J 



the artificial records, can be related to (AA/A) 2 
through the parameter a which is a known function 
of (AA/A) 2 [Lawrence and Jespersen, 1959]. 

For values of <r a as great as 3 radians, the linear 
relationship 

a d = 0.55 a a 

appears to hold. Figure 4 shows the observed 
relationship between (AA/A) 2 and P n . 



The development of these digital methods was 
sponsored by the Ballistic Missile Division of the 
United States Air Force, under D. 0. AF 04(647)-134. 
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